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ABSTRACT 

Models for the synchrotron emission of gamma-ray burst afterglows suggest that the magnetic field is generated 
in the shock wave that forms as relativistic ejecta plow through the circum-burst medium. Transverse Weibel 
instability efficiently generates magnetic fields near equipartition with the post-shock energy density. The 
detailed saturated state of the instability, as seen in particle-in-cell simulations, consists of magnetically self- 
pinched current filaments. The filaments are parallel to the direction of propagation of the shock and are about 
a plasma skin depth in radius, forming a quasi-two-dimensional structure. We use a rudimentary analytical 
model to argue that the Weibel filaments are unstable to a kink-like mode, which destroys their quasi-two- 
dimensional structure. For wavelengths longer than than the skin depth, the instability grows at the rate equal 
to the speed of light divided by the wavelength. We calculate the transport of collisionless test particles in the 
filaments experiencing the instability and show that the particles diffuse in energy. This diffusion marks the 
beginning of thermalization in the shock transition layer, and causes initial magnetic field decay as particles 
escape from the filaments. We discuss the implications of these results for the structure of the shock and the 
polarization of the afterglow. 

Subject headings: gamma-rays: bursts — MHD — instabilities — magnetic fields — plasmas — shock waves 



1. INTRODUCTION 

Gamma-ray bursts (GRB) afterglows have been ascribed to 
synchrotron emission from relativistic shocks in an electron- 
proton plasma. Detailed studies of GRB spectra and light 
curves have shown that the magnetic field strength in the 
shocked plasma (the downstream) is a fraction of eg ~ 
10~ 2 - 10~ 3 of the internal energy, while the energy in the 
emitting electrons is a fraction of e e ~ IP" 1 of the i nternal 
energy ( Panaitescu & Kumar 2002; Yost et al]|2003). Re- 
cently, Eichler & Waxman (2005) have shown that if only 
a small fraction / of the electrons are accelerated into the 
high-energy nonthermal tail, the observations can be fit- 
ted with values of e# and e e that are smaller by a fac- 
tor of / than the above estimates, as long as / is larger 
than the electron to proton mass ratio m e /m p . This still 
implies a magnetic field with eg > 10~ 6 . Furthermore, 
the measur ement of linear p olarizati on at the level of a 
few percen t ( Biornsson. Gud mundsson. & J6hannessonl2004t 
ICovino et al J 120041 and references therein') implies that the 
magnetic field in the synchrotron emitting region must deviate 
from isotropy. 

Simple compressional hydrodynamic amplification of a 
pre-existing magnetic field of the unshocked plasma (the up- 
stream) results in eg ~ 10~ 9 (Gruzinov 2001). Thus the req- 
uisite magnetic field must be generated in the shock itself 
or in the downstream. A leading candidate mechanism that 
produces a magnetic field near equipartition in th e shock 
transition layer is the transverse Wei bel instability {W eibel 
1959; Fried 1959), as was suggested by Gruzinov & Waxman 
i 19991) and iMedvedev & Loebl (I1999I) . 1 This instability and 
the magnetic field it produces are expected to play a crucial 
role in the thermalization of the upstream and in the shock 
dynamics. 

Recently, a great progress toward understanding unmag- 



netized collisionless shocks has been made by means of 
two- and three-dimens ional p article-in-cel l (PIC ) simulations 
(e.g. iLee&Lampd 119731 iGruzinovl l200lt ISilvaetalJ 



1 See Bret, Firpo, & Deutsch 1 2005) for a discussion of different modes of 
the Weibel instability. 



200a iFrederiksen et alJ ko04 iJaroschek et all l2004t 
Medv edev et alJ 120051 iKatol 120051) . Although the simu- 
lations do not resolve the e-p shock, they do show that within 
a layer ^100 proton skin depths wide, where the upstream 
and the downstream plasma interpenetrate, the transverse 
Weibel instability saturates. 

The saturated state of the transverse Weibel instability con- 
sists of magneti cally self-pinched current filaments (see, e.g., 
Fig. 1 in Frederiksen et al. 2004). The filaments are initially 
about a proton plasma skin depth in diameter and are paral- 
lel to the direction of shock propagation. As such they come 
with a magnetic field close to equipartition (eg ~ 0.1) that lies 
in the plane perpendicular to the direction of shock propaga- 
tion. However, the observed emission from GRB afterglows 
is expected to be produced at a dist ance of ove r ~ 10 9 plasma 
skin depths from the shock (e.g., lPiranll2005l and references 
therein). Thus, even if the filamentary picture correctly de- 
scribes the transition layer of GRB shocks, only the late evo- 
lu tion of these filam ents is relevant for the observed emission. 

Gruzinov (2001) pointed out that there is no obvious theo- 
retical justification for the perseverance of the magnetic field. 
Since the field forms only on a small scale — the plasma skin 
depth — one would expect that it also decays over the small 
distance comparable to the skin depth. On the other hand, 
based on a quasi-two - dimen sional picture of current fila- 
ments, Medved ev et alJ (|2005) suggest that the interaction be- 
tween neighboring filaments may result in a magnetic field of 
an ever growing coherence length, whereby e# is saturated at 
a finite value many plasma skin depths in the downstream. 

There are reasons to believe, however, that the quasi- 
two-dimensional picture is short-lived. For example, 
shock compression cannot be achieved in the region 
where the filaments make such an o rdered structure 
(Milosavlievic. Nakar. & Snit kovskvl 120 05). Lacking de- 
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tailed PIC simulations, however, the filament evolution, the 
isotropization of the upstream particles, and the production 
of shock jump conditions remain poorly understood. These 
processes, of course, are of immense importance for the mi- 
crophysics of GRBs and collisionless shocks in general. 

The two main processes that dictate the evolution the or- 
dered outcome of the transverse Weibel instability are the 
magnetic self-confinement of each filament and the interac- 
tions between neighboring filaments. Here, we use simple an- 
alytical arguments to explore the first one, namely the stability 
of a single filament. In §|2] we present an MHD model of such 
current filament. In §|3] we show that the filament is unstable 
to a kink-like mode, and estimate the growth rate of the mode. 
This instability destroys the quasi-two-dimensional geometry 
of the filament and produces a significant magnetic field in 
the direction parallel to shock propagation. In § |4] we study 
the motion of collisionless test particles in the background of 
a current filament undergoing the instability. We show that 
charged particles confined within the filament diffuse in en- 
ergy space. We suggest that this diffusion isotropisizes and 
thermalizes the upstream particles. In §|6j we discuss impli- 
cations of these processes for the magnetic field decay and the 
physics of GRB afterglows. 

2. EQUILIBRIUM WEIBEL FILAMENT 

Consider a magnetostatic equilibrium composed of an infi- 
nite cylindrical current filament extending in the z direction. 
Let r = be the center of the filament, where we use cylin- 
drical coordinates (r,9,z)- Let the current flow in the +z di- 
rection along the central axis. Further, assume that at some 
radius r = R, the current drops to zero. We refer to R as the 
radius of the filament. At r > R, there are currents flowing in 
other filaments in both directions (±z), so that the total current 
flowing through any given plane z = const vanishes. Since we 
are interested in the stability of a single filament, we ignore 
the effect on the filament of currents external to the filament. 

Particles giving rise to the current must move coherently, 
e.g., if there is a pinching azimuthal magnetic field Bg > 0, 
positive charges move in the positive z-direction, and nega- 
tive charges move in the negative z-direction. This is possible 
if the particles do n ot execute a full gyration in the magnetic 
field jAlfvenl!939h : fr actional gyra tions give rise to directed 
current as explained in Spitzer| i 19651) . This requirement im- 
plies that a filament's radial structure depends on its radius. A 
detaile d dis cussion of this structure can be found in lDavidsonl 
JT971 and iHondal ([2000); we here summarize the relevant 
aspects. 

Consider first the regime in which the radius of the filament 
R is smaller than about the plasma skin depth, R < 8. The 
skin depth equals S = C7 1 ' 2 /w p , where cj p = (47re 2 «/m) 1//2 is 
the plasma frequency, 7 is the kinetic Lorentz factor the parti- 
cles, m is the mass of the particles, and n is their density (in 
shocks, m is the electron mass; in e~p shocks, it is the proton 
mass). This directly implies that the current flowing through 
the filament / ~ ■KR 2 ne(3\\c, where j3nc is the average axial ve- 
locity of particles in the filament, does not exceed the Alfven 
critical current /a = -f f3\\mc 3 / e . Since the magnetic field is 
related to the current via Bg ~ I/Rc, the condition / < I a im- 
plies that the Larmor radius of the particle = ^f3mc 2 jeBg 
must exceed the radius of the filament r^>(j3 / ' (3\\)R > R- As 
a result, when R<8, all particles magnetically confined to the 
filament move in directed fashion. Here f3n depends weakly 
on r and thus the current profile across the filament is approx- 



imately homogeneous. 

However when R > 8, if the current were uniform within 
r < R and the average axial velocity of the current carrying 
charges were relativistic (/3m ~ 1), the Larmor radius would 
be smaller than the radius of the filament, and directed mo- 
tion of charges within the filament would be compromised. 
Directed current flow can still be maintained when R > 8 if 
the current is confined within a thin annular cylindrical region 
of width AR not exceeding the Larmor radius. The current 
then gives rise to a thin magnetic "wall" against which a par- 
ticle confined to the filament can be reflected. Then the Lar- 
mor radius within the current carrying layer can be written 
tl = 8 2 /AR, and the condition > AR implies that the cur- 
rent carrying layer is thinner than the skin depth, AR < 8, as 
seen in PIC simulations. In this regime, the saturated phase 
of the transverse Weibel instability consists current-carrying 
domains separated by thin magnetic walls. 

An alytic considerations ( Milosa vlievic. Nakar. & Spitkovskvl 
2005) and numerical simulations indicate that once the trans- 
verse Weibel instability saturates in collisionless shocks, the 
filament size is comparable to the skin depth. Motivated by 
these results we focus here on the case in which R ~ S. 

We assume that MHD equations apply and that the pressure 
tensor of the particles is isotropic. Both of these assumptions 
are oversimplifications. The pressure tensor near the axis of 
the filament is approximately isotropic when ~ R ~ 8. The 
purpose of the MHD model is to elucidate the physical mech- 
anisms and motivate specific collisionless PIC simulations of 
the saturated state of the Weibel instability. The simulations 
are the best way to test the theory over a range of parameter 
values. 

The toroidal magnetic field B = B(r)8 is related to the axial 
current density J = J{r)z via 

1 d 

4nJ= (rB). (1) 

r dr 

Within the MHD approximation, the fluid pressure P(r) satis- 
fies the equation of pressure equilibrium 

VP=JxB 

= -^^-(r 2 B 2 )r. (2) 

07rr z dr 

To construct an magnetostatic equilibrium filament, one can 
choose the radial dependence of the magnetic field, and then 
evaluate the current density and the pressure using equations 
Qand0. 

3. STABILITY 

It immediately follows from equation (j2ji that the fluid pres- 
sure inside the filament is larger than outside, the magnetic 
pressure accounting for the difference. This pressure imbal- 
ance is the origin of the unstable behavior that we explore 
below. The filaments are unstable to the well-known sausa ge, 
kink, and related MHD modes (see., e.g.. lHasegawal!9 75 and 
references therein) which result in the distortion of the fila- 
ment boundary. We here focus on a particular mode, the heli- 
cal kink instability, but expect similar stability criteria, growth 
rates, and particle transport in other related modes. 

Consider linear magnetostatic perturbations around the 
equilibrium described in § |2] Ac cording to the energy princi- 
ple, elegantly proven in Kulsrud ( 2005), a perturbation given 
by the Lagrangian displacement £ grows if the associated po- 
tential energy change is negative. The potential energy change 
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is given by ( B ernstein et all 19581) 

Q 2 



4tt 



+ J • (£ x Q) + 7f( V • + (£ • VP)( V • 



dx 

(3) 

where Q = V x (£ x B) is the magnetic field perturbation, and 
the volume integral extends over all space in the radial direc- 
tion and an averagin g over z-dir ection is assumed. 

Following Newcomb ( 1960), we consider a perturbation 
with displacement £ = (£ n £g, £.) of the form 



i(m6+kz) 



(4) 



where £,-(r), £e(r), and £ z (r) are real functions of radius. New- 
comb shows that non-axisymmetric perturbations (m / 0) that 
minimize W are incompressible, V • £ = 0, and have form 



6=- 

m 



2, A 



k 2 r 



d f i 



d 

dr k 2 r 2 + m 2 dr \ r 



ikr 3 



k 2 r 2 + m 2 dr 



(5) 



Note that the function £, identical to the radial displacement 
£ r , has been left unspecified. 

For perturbations of the form in equation Q, the potential 
energy perturbation in equation (0 can be expressed in terms 
of the radial displacement only, 



where 



7T 

w = - 



rm 2 Bl 



+ 8 e 



dr. 



(6) 



k 2 r 2 + m 2 
1 m 2 B 2 
r k 2 r 2 + m 



+m — 
dr 



B 



-m 2 B 2 e - 
r 

2 



B s d 
-2^-(rB e ) 
r dr 



(7) 



as shown in equations (16-18) of iNewcombl i 1 960ft (note that 
in our case B z = 0). 

The perturbations with m = 1 are special in that they need 
not vanish at r = and do not incur substantial cost in mag- 
netic energy as the field lines are bent only minimally. There- 
fore one is allowed to assume that the radial displacement is 
independent of radius, d^ /dr = 0, over the region with the 
non-vanishing magnetic field, which greatly simplifies the 
analysis. The true fastest growing mode may not have con- 
stant £, but one expects the constant £ approximation to come 
close. For these perturbations the radial displacements are 



£ r = £ = const, 

k 2 r 2 



& = i 1 + 



k 2 r 2 +l 



fa k 2 r 2 + V' 



(8) 



With this, the energy per unit length along the z-direction be- 
comes 



Equation l|9} tells us that all wave numbers k > are unstable. 
In the long and the short wavelength limits, the energy can 
approximately be written as 



W 



-\-Ki 2 k 2 R 2 B 2 , (kR^il). 



(&?>1) 



(10) 



where Bg is the RMS magnetic field, and Be is the magnetic 
RMS field averaged per unit log-radius. 

Equation (|9j tells us that all m = 1 modes with finite wave- 
lengths (k > 0) are unstable. An estimate of the linear growth 
rate is given by 

r~w 

r~J--, (11) 



K' 



where 



1 



K = - I tt&+$ + g)Td6dr. 



(12) 



The value of K for perturbations in equations ^ equals 



K = 



[5k 2 R 2 -3ln(k z R z +Y)] 



2 D 2 



4k 2 



(13) 



where p is the average mass density. In the long and the short 
wavelength limit, therefore, 



r 



kBe/yip, (JfcR « 1) • 



(14) 



Note that the growth rates are proportional to the "nonrela- 
tivistic" Alfven velocity va = Be/ \/4np multiplied by k and 
ST 1 , respectively, in the long and the short wavelength limit. 
This is expected since the instability is driven by a pressure 
imbalance. The growth rate in the long wavelength limit can 
also be expressed as T ~ (7r/2) 1 / 2 /3||^w p , where as before uj v 
is the plasma frequency. 

Note that the magnetic field acquires a component parallel 
to the axis of the filament. The RMS strength of the parallel 
field is B z = Q z r ^ \k^Be in the long wavelength limit. 

The above analysis is valid only as long as the perturbation 
is nonrelativistic (i.e., £ < c/T) and the wavelength of the per- 
turbation is larger than the Larmor radius. Since we here con- 
sider the filaments with r^ ~ R ~ 5, the analysis is valid only 
for wavelengths with 2ir/k > r L > 5. In such filaments the 
growth rate is T ~ (3\\ck. 



4. TRANSPORT IN AN UNSTABLE FILAMENT 

We next address the orbits of collisionless test particles in 
the electromagnetic field defined by the magnetostatic equi- 
librium (§0 and the linear unstable perturbation (§0. Con- 
sider a particle with unperturbed orbit confined to the region 
r <R, and assume that the motion is directed along the axis, 
dz/dt > 0. The momentum of the particle p is then a periodic 
function of z, as are its radial excursion and azimuth. The 
magnetostatic equilibrium is electrically neutral and the en- 
ergy of the particle, E = ync 2 , where 7 = (1 +p 2 /m 2 c 2 )^ 2 , 
subject to the magnetic field of the equilibrium only, is a con- 
stant of motion. The time variation of the magnetic field of the 
unstable perturbation induces an electric field given by Ohm's 
law 

VxB 

E = , (15) 

c 

where B is the magnetic field of the equilibrium, and V = 
d£/dt is the bulk velocity associated with the perturbation. 
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In long wavelength perturbations, the displacement is mainly 
perpendicular to the z-axis and the electric field is thus mainly 
parallel to the axis. 

In writing equation (I15> . we have assumed infinite conduc- 
tivity, and have ignored terms of the form J x B, dj/dt, and 
VP. These terms can be neglected to the first order in k5 and 
T / uip, i.e., for wavelengths longer than the plasma skin depth. 

Kinetic energy of a test particle with charge q incurs a 
change of 



A£ = q [ ^dt' 



7m 



(16) 



To the first order in the perturbation, the energy change equals 



A£=-qT / /?•(£ xB)rfr' 

J — OO 

= -?r / [&A.sin(jfcz+0) 



+£ r /? z cos(fc:+0)]e Bgdt' 



(17) 



where B = p/jmc and as before T denotes the growth rate of 
the instability. 

The energy change can be evaluated for the incompressible 
perturbations given in equations (|8j. The displacement of the 
fluid is an oscillatory function of the monotonic variable z(f), 
as is the displacement of the particle. Resonances occur when 
B\\ck ~ uig, where Bnc is the average axial velocity of the par- 
ticle, and cue is its azimuthal angular frequency, which will 
differ from one particle to another. At a resonance, the parti- 
cle experiences a coherent electric field over most of its orbit. 

5. THERMALIZATION 

To explore the effect of energy change derived in equa- 
tion 1171 . we consider the simplest Maxwell- Vlasov equi- 
librium of a magnetically self-pinched cu rrent filament 
jHammer & Rostokerl 1 197(1 iDavidsoi] 1 1974ft . in which all 
particles have the same total energy £ = ^mc 2 and the same 
axial canonical momentum V z = (3^ mc, both of which are 
constants of motion. Here, 70 is the Lorentz factor of the par- 
ticles and < B11 < (l-% 2 ) ] ^ 2 is a constant. The canonical 
momentum equals V = p- (e/c)A, where A is the vector po- 
tential. 

Particle density n inside the filament associated with this 
equilibrium is uniform. The equilibrium has phase-space dis- 
tribution function 



K£,T>zY- 



27T7om 



S(£- jotnc 2 )5(P z - B\\ jomc), 



(18) 



where S(x) is the Dirac delta-function, not to be confused with 
the skin depth. The magnetic field and the vector potential 
inside the filament are given by 



B„ = 



A,: 



/3||7 rac 2 

7 
/3||7 mc 2 



Mi 



(19) 



where as before S is the plasma skin depth inside the filament, 
and I„(x) is the modified Bessel function of the first kind. The 
filament radius R is the solution of the equation 



1 



7o 2 /?f 



(20) 




.01 



kR/27T 



FIG. 1. — The RMS lineal' energy change A£ of a particle in units of 
(£/R)"/mc 2 as a function of the wave number k for the model described in 
§g]with /3|| = (0.05,0.5,0.95) from top to bottom and 70 = 10. Here, £ is 
the final value of the displacement associated with the instability discussed in 
ij l3l The energy change was calculated using Monte-Carlo integration. 



Although extremely simple, this model correctly captures 
many of the characteristics of current filaments seen in PIC 
simulations, including the concentration of the current flow 
near the edge of the filament when R > 5. 

Figure [2 shows the RMS linear energy change of parti- 
cles drawn from the distribution of equation (II 81 with Bu = 
(0.05,0.5,0.95) as a function of the wave number k of the 
unstable perturbation. The magnetic field of the equilib- 
rium filament is given by equations ( 1191 . and the filament 
is undergoing instability with displacement given in equa- 
tions (|8} and growth rate defined by equations l|9), M Q , and 
\\2\ . For Bu > 0.5, the linear energy change is maximum for 
kR ~ (0.1-0.4) x 2tt. It should be kept in mind, however, 
that the pressure tensor in filaments with Bu <C 0.5 and those 
with 1 — B\ <C 0.5 is anisotropic, and thus the growth rates 
calculated in §|3]may not be accurate. 

The instability affects in a similar way the pitch angle of the 
particles; the pitch angle can be defined as the angle subtended 
to the central axis at the instance of closest approach to the 
axis. We have focused on the energy because it is a simpler 
calculation. 

As the orbits of particles inside a filament are perturbed by 
the instability, some of the particles that are initially confined 
to the filament escape. For example, consider a particle mov- 
ing in the meridional plane (9 = const) that has radial velocity 
B r when crossing the central axis. Assume that the particle 
is marginally confined to the filament, i.e., B 2 ~ R/rt, where 
tl ~ jmc 2 /eBg is the approximate Larmor radius. Evidently, 
an increase of the pitch angle tfi = sin _1 [/? r (l -7~ 2 ) -1 ^ 2 ] at a 
constant 7, or an increase of the energy £ = -fine 2 at a con- 
stant pitch angle, both caused by the perturbed electromag- 
netic field of the instability discussed above, can liberate the 
particle from the filament. 

A filament experiencing a growing kink instability contin- 
uously sheds particles in this fashion. The liberated particles 
have sufficient energy to visit neighboring filaments of either 
sign of the current, and can be thought of having joined a 
thermalized pool. The loss of particles implies decrease of 
the current flowing through the filament, which in turn im- 
plies the decay of the magnetic field. Only the particles that 
remain confined to the filaments contribute to the coherent 
toroidal magnetic field of the filament. Rapid field decay can 
be partially offset by an increase in the current per particle by 
electric fields induced during flux loss; robustness of the cur- 
rents depends on the competition between the scattering and 
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the induction. 

6. DISCUSSION 

We have argued that the quasi-two-dimensional structure 
of the transition layer in collisionless shocks is disrupted by 
pressure-driven instabilities. We have also shown that the dis- 
ruption is accompanied by a redistribution of particle ener- 
gies, which can be interpreted as the onset of thermalization 
and magnetic field decay. Long term evolution of the mag- 
netic field must therefore be addressed in the context of the 
three-dimensional turbulence that ensues after the Weibel fil- 
aments have been disrupted. The latter regime remains poorly 
understood. 

An alternate interpretation of the early decay of the mag- 
netic field refers to the hier archical merging of current fila- 
ments due to Lor entz forces (Gruzinov 2001 ; Med vedev et al.l 
200.4 iK^tol 12005). It can be asked whether the merging or 
the pressure-driven instabilities prevail. Since the maximum 
current that can flow through a filament is limited (see §|2j, 
magnetic energy density must eventually decay under merg- 
ing. However, PIC simulations of <? ± shocks (Spitkovsky & 
Arons 2005, private communication) show that the current in- 
side a filament is tightly shielded by a reverse current flowing 
just outside the filament, i.e., opposite current filaments are 
tightly packed and most of the current flows near the edge. 
The shielding currents reduce the Lorentz attraction and slow 
the growth of the magnetic field correlation length via merg- 
ing. 

Meanwhile, the filaments are susceptible to the instabili- 
ties independently of the shielding. Therefore, we expect the 
dynamics of single filaments to be governed by the instabili- 
ties. In particular, the instabilities may drive merging between 
the filament fragments in three dimensions. The highest- 
resolution published simulatio ns of cold shell co llisions with 
two particle species by iFrederiksen et al.l (|2p04) show clear 
evidence for progressive bending and kinking of the "proton" 
(m p /m e = 16) filaments (see their Fig. 2). 

Similar mechanism was recently studied by 
IZenitani & Hoshinol (120051) . who carried out two-dimensional 
PIC simulations of the relativistic drift-kink instability in an 
infinite e current sheet confined between reversed magnetic 
fields. As the sheet bends in the direction perpendicular to the 
magnetic field, an alternating electric field is induced parallel 
to the sheet. The authors compare the growth rates measured 
in the simulations with predictions from two-fluid theory 
and find good agreement for kX < 0.7, where k is the wave 
number and A is the thickness of the current sheet. They also 
find that in the nonlinear stage of the instability, the electric 
field becomes coherent in the central region of the current 
sheet, which accelerates particles and efficiently dissipates 
the magnetic energy energy. 

If external shocks in GRBs resemble the observed outcome 
of relativistic shell collisions in PIC simulations, the instabil- 
ity discussed here has implications for the interpretation of 
the observed weak linear polarization of the afterglow. The 
optical emission of GRB afterglows is linearly polarized at 
a level of a few percent, implying that the magnetic field 



in the emitting region is anisotropic. There are two differ- 
ent forms of magnetic field anisotropy that can produce this 
polarization. The first is a field parallel to the plane of the 
shock that is coherent on scales exce eding the plasma skin 
depth by many ord ers of magnitude (G ruzinov & Waxmanl 
1999; Granot 2003). The second is a combination of a ran- 
dom magnetic field within the plane of the sh ock and a 
non-a xisymmetric geometry of the emitting region (|Gruzinov | 
[T9^|Sar]|T999l:lGhisellini & Lazzat|1999llGranot & Konigll 
2003t lNakar. Piran. & Waxmanl2003tlRossi et all2004l) . 

A coherent magnetic field within the plane of the shock is 
unlikely to be produced in the shock itself since there is no 
preferred direction within the plane of the shock. Such field 
could in principle result by amplification of a pre-existing, 
ambient coherent magnetic field. This possibility, however, is 
implausible because the field in the medium into which the 
GRB ejecta plow is expected to be far too weak. The sec- 
ond possibility is that the coherence length of the magnetic 
field generated the shock grows rapidly. Since at any given 
time the observer sees many causally disconnec ted regions 
JGruzinov & Waxmanl 1999tlNakar & Orenl2004l) . the coher- 
ence length of the downstream magnetic field should grow at 
a rate close to the speed of light in order to produce a polariza- 
tion at the level of a few percent. We do not know of a mech- 
anism that would facilitate such growth, especially given that 
the Alfven speed in the weakly magnetized relativistic down- 
stream plasma is much smaller than the speed of light. 

An anisotropic field, but random within the plane of the 
shock, appears to be a much more appealing possibility be- 
cause the shock breaks the isotropy and could in principle re- 
sult in a different mean field strength in the parallel and the 
perpendicular direction. In this case the level of polarization 
also depends on the degree to which the axisymmetry of the 
emitting region is broken. Assuming that the magnetic field 
lies entirely in the plane of the shock (i.e., B z = 0), and that 
the geometry of the emitting region is as expected during the 
jet-break in the light cur ve, the level of p olarization is at most 
20% -30% dSarilll999t [Rossi et all 120041) and in some sce- 
narios it can be as small as a few percent. If B z « Bxy, the 
polarization level is significantly reduced. Therefore, the typ- 
ical observed polarization of a few percent requires that either 
B, < B„./2 or B z > 2B xy . Even this minor difference between 
B z and B xy presents a theoretical challenge in view of our re- 
sults, which suggest that B z arises quickly within the transition 
layer as the current filaments associated with the small-scale 
magnetic fields are effectively destroyed. We speculate that 
the magnetic field quickly becomes isotropic in the rest frame 
of the shocked plasma. 
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